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Abstract 



The averages of ratios of characteristic polynomials det(A — X) of N x N random 
matrices X, are investigated in the large N limit for the GUE, GOE and GSE ensem- 
ble. The density of states and the two-point correlation function are derived from 
these ratios. The method relies on an extension of the Harish-Chandra-Itzykson- 
Zuber integrals to the GOE ensemble and to supergroups , which are explicitly eval- 
uated as solutions of heat kernel differential equations. An external matrix source, 
linearly coupled to the random matrices, may also be added to the Gaussian distri- 
bution, and allows for a discussion of universality of the GOE results in the large N 
limit. 
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1 Introduction 



In this article we consider correlation functions involving characteristic polynomials 
of N X N random matrices X belonging to the Gaussian unitary (GUE), Gaussian 
orthogonal ensemble (GOE) or Gaussian symplectic ensemble (GSE), of the following 
type 

F(A„---,A,;/.,---/..) = (n f^|^"~^| ). (1) 

a=i det(/i„ - X) 

The reason for considering expectation values of such characteristic polynomials are 

the following. First they turn out to be simpler than the usual correlation functions 
k ^ 

of the type (J]^ Tr — ), but they may be used to recover the same information. 

For instance 

d . det(A-X) ^ . (0) 

and so on. Next we know ^ that correlations of products < na=idet(Aa — 

X) > are universal in the Dyson limit, in which the distance between the A's is 
of order of the mean spacing |jl], ^ ^. This universality is known to hold if the 
probability distribution is non-Gaussian, or if it includes an external matrix source 
linearly coupled to the random matrices. Presumably this extends to the ratios (0) 
considered here and some evidence in that direction will be presented below, when 
the Gaussian measure if modified by the presence of an external source. 

The calculations will be shown to involve integrals over groups which go often 
beyond the case considered by Harish-Chandra, Itzykson and Zuber (HIZ) [^, |10[and 



the main point of this article is to analyze such cases. Such integrals appear at 
many places in random matrix theory. For instance they are essential in the problem 
of a random Hamiltonian H, an Hermitian N x N matrix, which is the sum of a 
non-random Hq and of a Gaussian random V : 

H = Ho + V (3) 

with 

PiV) = le-fT^^^ (4) 



The probability law for the matrix elements of H is thus 

When Hq is non-zero this measure is not invariant under a change of basis, i.e. under 
the orthogonal, unitary or symplectic groups appropriate to real symmetric, hermitian 
or quaternionic matrices. If one wants to write now the probability distribution for 
the eigenvalues of if, one writes H = UAU~^ with A diagonal, and U a group element 
appropriate to the ensemble there is a well-known Jacobian 

N 

dH = dU\A{X,,---,X^)fl[dXa, (6) 

1 

with /? = 1, 2 or 4. and the probability distribution for the eigenvalues is given by 
P{Xu ■■■,Xn) = || A(Ai, ■ ■ ■ , A;v) fe-TE^'^ J rf[/e^^(^^^-^^o). (7) 

Therefore a group integration has to be performed. 

Similar integrals occur in many places of random matrix theory. For instance for 
chains of random matrices on a lattice with couplings of the form expTiHnHn+i, 
clearly the joint probability for eigenvalues involves a similar group integration over 
the relative U~^Un+i- 

In the present work the correlation functions (|l]) will be written as supersymmetric 
integrals and diagonalization through supergroup transformations will lead us to 
consider new types of HIZ group integrals. 

The literature on supersymmetric integration is very rich. Efetov has already 
obtained expressions of the correlations for those same ratios ([T|) by supersymmet- 
ric techniques. In the unitary case, in which X is a complex hermitian matrix, the 
calculation of the n-point correlation is done explicitely with the help of the HIZ inte- 
gral, which has been extended also to the supermatrix case |TT|, . However for real 
symmetrix matrices, for general k, the supermatrix formulation becomes complicated 
since the relevant HIZ integrals are no longer explicit and simple. Recently, Guhr 
and Kohler |]13|, H have studied those integrals by the recursion formula of the size 



N of the matrix, but it is a different approach which is proposed below. In a previous 
work, we have considered such integrals for the symplectic group , with the help of 
differential equations satisfied by the heat kernel for motions in the orthogonal group. 
This approach may be extended to include an external matrix source, a modification 
for which the standard orthogonal polynomial methods do not apply any more. 

Since we have considered such HIZ integrals in the past, let us summarize what 
was known earlier, beyond the standard HIZ integral for the GUE. In the GOE case a 
mapping of the N x N integral for (Hi det(AQ, — X)) to an integral over kxk matrices 
invariant over the symplectic group was first derived. Then one found, through the 
heat kernel differential equation, that the HIZ integral was exactly given in that 
case by a semi-classical approximation corrected by a finite number of terms. Those 
additional terms are important, since they are of order one in the Dyson limit ; once 
those terms are known one can determine this limit by a saddle-point method. Here 
we return to those heat kernel equations for the supergroups relevant to our present 
problem. It turns out that now the series of correction terms to the semi-classical 
limit does not terminate. However in the scaling limit this series may be determined 
explicitly, and the problem at hand is thus solved. 

The set-up of this article is as follows : We first consider the average of a single 

ratio (- — [ \) for the GUE and reduce it directly to quadratures, or alterna- 

det(/i — X) 

tively we make use of supersymmetric techniques with a simple version of the HIZ 
integral. Then the same is done for the GOE. For a single ratio the integration 
over the supergroup variables is done through a generalization of the HIZ formula, 
which is exactly given by the semi-classical approximation plus a finite number of 
corrections. We then apply those supersymmetric techniques to higher correlation 
functions. The required HIZ formula for supergroup integration leads to a differ- 
ential equation whose solution shows than the semi-classical integration has now to 
be corrected by an infinite series. This infinite series originates from the part of the 
supergroup integration which comes from the orthogonal sub-group, and each term of 



3 



this series is of order one in the Dyson hmit. However in the hmit of large matrices it 
turns out that one needs only to consider a special case of those HIZ integrals, when 
the source matrices have only two distinct eigenvalues. Then the whole series may be 
found in closed form and one obtains explicitely the large N limit of the correlations 
functions for ratios of characteristic polynomials. This is then generalized to k-point 
functions and, continuing in k, one can study the zero-replica limit /c — 0, which is 
an alternative way of deriving the usual correlation functions from those determinen- 
tal expectations. The Gaussian symplectic ensemble (GSE) is discussed within the 
same appraoch. Those techniques are then extended further to probability measures 
involving an external matrix source. This allows one to verify that the results in the 
Dyson limit are not modified by this source. Level spacing distributions are then 
derived within the same methods. The case of arbitrary (3 is then examined. Finally 
edge singularities are discussed as well. 

2 Single ratio of characteristic polynomials in the 
GUE ensemble 

Let us first consider a single ratio of characteristic polynomials, defined as 

The averages < ■ • ■ > are computed with respect to the Gaussian distribution 

P{X) = ^expi-^TrX') (9) 

where the matrix X is a x complex hermitian matrix. 

The determinant det(A — X) may be expressed as an integral over Grassmann 
variables 6 and 6 0. The determinant det(;U — X) in the denominator, may be 
expressed as a Gaussian integral over commuting variables z and z*. Then 

TV 

FAr(A,/i)= / '[[dzldzad9ad9a<expiN[9a{XSab-Xab)9b + zl{fi6ab~Xab)zb]> (10) 

a=l 



in which /j, is complex, with a small positive imaginary part. The normalization for 
Grassmann integration which is used here is 

/ deS 99 = - (11) 

We also adopt the convention 



(a/3) = ap (12) 

which implies 

a = —a (13) 

in order to maintain 99 invariant under the bar operation. 

Then the Gaussian integration over X is easily performed since 

with Yji = —9i9j — z*Zj. This gives 

TrY^ = TrYijYji 

= -{9A)^ + {z*Zi)'' + 2{9iZi){9jz*) (15) 




TT 

N 



One then introduces auxiliary commuting as well as Grassmann variables 

j ^-Nfjv+Nr,{ez)+Nfj{z*e)^^^;^ ^ ^-N{ez){z*e) (ig) 

This allows to express Fjv as 

X [ J dzdz*d^d^e^(^^-^')^"^+^^('^-'')^*^e^^^^*+^''^"^ ]^ (17) 

Two strategies are possible at that stage : either keeping both the commuting and 
Grassmannian degrees of freedom, and use the supersymmetry of the integral, or 
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integrate out explicitely the Grassmannian variables. They are equally simple for 
the single ratio of characteristic polynomials considered here, however for higher 
correlation functions the use of the invariance under supergroup transformations turn 
out to be more powerful, 
(i) integration over the Grassmann variables 

Let us first give the expressions that one finds if one first perform the integrations 
over 6, 6, z*, z. This yields 



2 J u — 



rjT] 



jj, — b {jj, — by 

^2n'J^ifi-b)^ (^_6)iv+iJe ^^^^ 
It is convenient to use the (appropriately normalized) Hermite polynomials 



^ / (A + it)"e-* dt (19) 



which satisfy the recursion formula 

Hn{\) - 2\Hn-i{X) + 2{n - l)if„_2(A) = (20) 
Then one obtains from (ITBI) 



F.(A,.)^^[i^.(A)/_ db---^-H...iX)J_ dt---^] (21) 



in which fi has a small positive imaginary part and we have used for convenience 

A = Ayy ^' = ^'\^■ (22) 

After repeated integration by parts this gives 

F^(A,/i) = ^^^^^_^^J _^ db^—^[HM{X)Hr,-,{b) - H^.^{X)H^{b)] (23) 
from which follows for an infinitesimal positive imaginary part of /i 



ImF;v(A,^) = - ^^^^'^^^, e-^'[g^(A)if^,i(/i) - H^^^iX)H^i^)] (24) 



which vanishes as expected for A = /i since F^ifi,^) = 1. 
The density of states p(A) is thus given by 

p(X) = -\im±-^lmF^iX,f^) (25) 

and one recovers the well-known finite expressions 
(ii) use of supergroup symmetry 

Let us return to (|1^ and use the supergroup structure. One first write 



- b'ee - ibz*z + f^Oz* + r^Oz = -iiz*,e) Q ( ^ ) (26) 



with 

n-f 

—irj —ib' 



Q=[ L Z) (27) 



Then 



istrg2 = b^ + ^ 2f]7] (28) 



in which we have used the supertrace notation "Str" defined as follows. If one de- 
composes the supermatrix Q into four blocks 

in which A and D have commuting matrix elements , whereas the matrix elements 
of C and B are anticommuting. Then by definition 

StrQ = trA - tiD. (30) 

Similarly one defines the superdeterminant 

det{A- BD~^C) detA 

The notation is justified by the fact that 

J dededz*dze-'^^^ = (Sdet[Q])-i 
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if(f)= (^^j and 0= {z*,e)). If one defines tlie diagonal supermatrix 

(32) 

This leads to the compact expression for Fn{X, fi) as an integral over the supermatrix 

Q, 



(Sdet(Q-A)) 



N 



One now performs a superunitary transformation to diagonalize Q 

ad) (34) 



9={1 !) (35) 



with 

,7 d 

in which (3 and 7 are anticommuting. Then, with the sign convention used here [z*, 9) 
transform by 

«'-(;* JO- p'') 

Unitary transformations g'^ g = 1 leave invariant the quadratic form z*z + 99. One 
may diagonalize the matrix Q by such a unitary transformation 

Q-o{l °)«* (37, 

Then the integration measure on the matrix Q may be replaced by an integra- 
tion over the eigenvalues ui,t and a group integration, up to a Jacobian which is 



proportional to — it) . The representation (33) requires then to consider the 



supergroup equivalent of the HIZ integral defined by 

^g^NStrgQg-^K (38) 



where A and Q are diagonal matrices. 



This integral yields 

N 



A)e^("^-'*^). (39) 



27r2 

The derivation of this formula relies on a formalism of differential equations for / 
which will be presented in detail in the coming sections. Note that the integral (^) 
vanishes when the matrices A or Q are multiple of the identity ; this is a consequence 
of the Grassmannian degrees of freedom which imply that J dg x 1 = 0. 
Therefore we have obtained the representation 

fD — -^(/^^ — '^^) POO POO /}-{- 1 

zh'' J-oo J-oo u u — it 
The density of states is then recovered through 

p(A) = -Imhm^ r-/-rft.(— )^^e-f(-^+*^)-^*^+^-'^ 



/i^A d^i J-oo 2-K J 2n u u ~ it 
1 dt r du -it 1 ^-Kiy?+t^)-im\+Nu\ ^i) 
N J-oo 2t: J 2iTi^ u ' u - it ^ ' 

where the contour integral of u circles around the origin. This exact representation, 



valid for any N , has been derived earlier by a completely different method |15 
without supersymmetry. The semi-circle law is recovered from there easily, in the 
large-N limit, by a saddle-point integration (it is simpler for that to consider dp/dX). 

3 Single ratio of characteristic polynomials in the 
GOE ensemble 

We now consider a single ratio of characteristic polynomials of GOE, defined as 

^ s det(A-X) , , 

The averages < ■ ■ ■ > are computed with respect to the Gaussian distribution 

P{X) = ^exp(-|TrX2) (43) 
where the matrix X is a x real symmetric matrix. 



As for the GUE, the ratio in ( ^21 ) is expressed by 



N 



Fjv(A,^) = / n dzldZadOadOa < e^^lN[ea{\5ah-Xab)0b + zl{^l5ab-Xab)Zb\ > (44) 
a=l 

in which /i is complex, with a small positive imaginary part. We use the normalization 
for Grassmannian integration of the previous section. 

Then the Gaussian integration over X is easily performed since 

with Yji = —9i9j — z*Zj. This gives 

TrY^ = TrYijYji 

= -m)' + {z*z,)' + 2{e,z,){e,z*) (46) 

TrYY^ = 2{e,z:){z^e,) + (z*zn(zjZj) (47) 

One then introduces auxiliary commuting as well as Grassmann variables 
IN r 



N 



TT J 



2N 

TT 



e 



-2Nn'v'+Nri'(ez'')+Nn'{ze)^-i^^, ^ ^-!L{ez'){ze) 



2N 

This allows to express as 

Fr,{X,fi) = ^Jdbdb'dudu*dr]dvdr]'dv'e-''^'"+'''"+^*^+^'^'^+^^''^'^ 
X [ J dzdz*ded9e''^'^-''^''-^'''^>'-'^''' 

Nf)ez'+N-qez+^N{u*z'^+uz*^)+Nri' z'e+Nf)'ezT^N 



(49) 



Two strategies are possible, as explained in the previous section. Let us first give the 
expressions that one finds if one first perform the integrations over ^, 9, z* , z. Noting 
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that 

the integration over 2;, 2;*, ^, 6*, yields 



((/i-6)2- |u|2)l/2 

{wf/'r/ + u*r7r7' — (/i — 6)(r7r7 + r^'V)}]^ (51) 



One can next integrate over the remaining Grassmann variables f],ri,fi' ,■)]' and 



denoting D = J (n — by — one obtains at the end 



X 



^-N{h'+b'^+u*u)^^^y^^^^* (52) 



It is convenient to use the (appropriately normalized) Hermite polynomials defined 
by (|T9|) . Then, using a few integrations by parts, one obtains from ( ^2]) 

1 _ f+CO (jh _ roo 1 

+ ^i/._.(A)/_^e-MA-2*)/ <ipe-|^_^j,_j^ (53) 
in which /x has a small positive imaginary part and we have used for convenience 

A = Av^ /2 = fiVN. (54) 



The relation ( pHD for the average ratio allows one to find in particular the density 
of states, and the density of states p(A) is thus 

p(A) = - hm ^ AiniF^(A,^) (55) 

when ImA 0. 

The explicit formula (|53|) allows one to recover standard results derived by the 
method of skew orthogonal polynomials 0. In order to write explicitely the imaginary 
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part of Fn one has to distinguish between even and odd N . Let us for instance work 
with N = 2M and check (||) for N=2 and N=4. Using 



+00 Q-b ^ 



00 



and 



Im / dp- — = -vr sgn(/i - b)e-^'^-''^ (57) 







where sgn(x) = +1 for x > and —1 for x < 0. One thus finds, for N = 2 

v^ImF2(A, /x) = (yu - A)e"^' - A(/i - A)5(/i) (58) 

in which 

B{x) = e-'''/^ r dye-y''\ (59) 

JO 

For = 4 we use , 

foo 1 1 poo 1 

e-^rfp = + / e-Pdp (60) 

and obtain 



{{p-hY-pf ' (^-6)2 Jo {ii-hY-p 



V^ImF4(A,/i) = 4e-^'[^2(A)if3(p)-i^3(A)if2(/i)] + ^(/i-A)if3(A)(2/ie-^' + i?(/i)) 

(61) 

We see that, as it should, ImF^^p, p) vanishes (since Fn{p,p) = 1). The above 
expressions agree with those derived from orthogonal polynomials; 



^2(A,/i) = 0o(A)0o(p) - 0o'(A) / Mt)dt 

Jo 

= 0o(A)0o(p) + 0i(A)0i(yu) + 0i(A) r02(t)c?t (62) 

Jo 

where 0„(A) = (2"n!v/^)"^/V^'/2^„(A). We find ImF2 = {p-X)S2{X, p)e-'''/^+^'/^ 
The S2{\, p) is the diagonal part of a kernel, which is a 2 x 2 quaternion matrix [0]. 
For general N, Sj\f{X,p) is given by 

^^(A,/i) = ^ 0,(A)0,(/i) + (-)i/20;v-i(A) / 0^(t)cit 



i=0 

AT 



E 02^(A)02.(P) - E <^'2.(A) r02.(t)rft (63) 
j=0 i=0 
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In the large N limit, the second term in (pSD in the first line is negligible compared 
to the first one. Indeed since 0Ar(t) is a rapidly oscillating function, its integral in 
the second term is asymptotically small. Then one recovers the semi-circle law in the 
limit A — > /i. Therefore the first term of 5'iv(A,/i) is dominating ; it coincides with 
the well-known universal kernel for the GUE ensemble 

Let us further check for N=4 case. In this case, we obtain 

lb 

+ [l^iH,iX)+^-iAf,'-l)H,{X)-^f,'H,iX)]e-'' (64) 

Then we have 

hm /^ImF4(A, /x) = -(/x^ - - (S/x^ + ^)e-^^ (65) 

The right hand of this equation agrees precisely with 5'4(/i, /i) in (^) . This allows 
one to check also the density of states since 

d 

lmhm—FN{X,fi) = -S{fi,fi) (66) 

We have verified the known result of the density of state in the GOE, p(A) = 
^^(A, A). 

We are interested in the large N limit and the choice of variables b, b' and p 
hereabove, is not the best one. Let us define instead ui, U2 and t, (analogous to the 



variables u and t which have been used for the unitary case[I^, |I6[ ; a comparison 
with the unitary case is very useful) as follows 

U1+U2 {ui - U2f 

b= — ^ ,P = ^ (67) 

and b' = t. The geometric meaning of those variables will become clear later when 
the formalism of graded supermatrices is introduced. Then we have from 



27r2 J J J-00 ' " (MiM2)^/' 



Ml — U2\ 



^ ^ 2it{uiU2) AN {it)^uiU2 
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This representation is well adapted to the large-N limit. The saddle points are 

-iX ± 

t" = = t+,t_ 



case 



Ml = = U+,U- 

= '^^'f^ . (69) 

We have analyzed before a similar large N behavior for the kernel of the unitary 
1% . The factor {it/^/u\U2)^ , is oscillatory in the large N limit unless we take 



the pair of saddle points with {t,Ui,U2) = (/!:+,«+,«+) and (t-,U-,U-). The other 
six possibilities are sub-leading. The fluctuations around the saddle-point give 

1 1 



(70) 



(/"(t.))V^ ^ 

and for Ui and U2 we have similar Gaussian fluctuations to take into account, together 
with the factor — M2I. Then the Gaussian integration for the u- fluctuations yields 



1 



The second term in the integrand of (|68| ) becomes at the saddle point 



(71) 



1 Ui + U2 1 1 1 1 , . 

1 -I L_! £ _| = 1 H I (72) 

2it{uiU2) 4:{ityuiU2 itcUc '^{itcYul 

which is (1 + 2it~[^^^^ Noting that Uc = itc, we find that this factor cancels the results 



(|70D and (^) of the Gaussian fluctuations. In the large N limit, then we have 



ImF/v(A, /i) = sin(A^(A - ^i)V2 - X^) (73) 



piX) = ^lhn|sm[iV(A-;.)v/23A^] 

= -V2^ (74) 

TT 



as expected. 
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4 Supermatrix diagonalization 



We have obtained an integral representation of the averaged ratio of characteristic 
polynomials (^) by a Gaussian integration over ^'s and z^s. If, instead of integrating 
out the anticommuting variables, we keep bosonic and fermionic (supersymmetric) 
degrees of freedom, the formulation becomes more transparent. Furthermore it turns 
out that the diagonalization of the supermatrices by a supergroup element provides 
an extremely useful representation in terms of the eigenvalues. 
Let us define the super-matrix Q by 



Q 



( ^ 

-u 
-ir] 
\ if]' 



—u 
h 

-irj' 
if] 



iri' \ 
ii] 

■ib' 
-lb' J 



IT] 

ifj' 



(75) 



which is such that one can write the various quantities which appear in the integral 
representation (|^) as follows : 



-ibz* z + -{u* z"^ + uz*'^) 



f]ez* + rjOz + ri'z*e + f]'ez - b'ee 



iz\z,e, 



Q 



e_ 



and 



1 



StrQ^ 



b^ + b'^ + u*u + 2r]T] + 2r]'r]'. 



From the definition (^) of Sdet[Q] , one has 

6^ - 

Sdet[Q] = 



u u 



(76) 



(77) 



(7^ 



[b' + (6(?7?7 + r]'r]') + urj'T] + u*riri')]'^ 

Similarly we define the diagonal supermatrix A , with elements (/x, /i, A,A) on the 
diagonal. This leads to a compact expression for F/v(A,yu) as an integral over the 
supermatrix Q, 



dQe 



'-StrQ'^-iNStrQA 



(SdetQ)^/2^ 



(79) 
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which is identical to (pBI). 

Superunitary transformations on the variables are defined by 







/ z'\ 


z* 




z'* 


e 


= 9 


6' 


\o ) 




\0' J 



with 



9 



a b 
c d 



in which h and c are 2x2 matrices with anticommuting elements. Then 



9 



c 

-b d) 



12) 



Unitary transformations g^g = 1 leave invariant the quadratic form z*z + 99. One 
may diagonalize the matrix Q by such a unitary transformation 



Q = 9 



U2 



it 



V 



9' 



(83) 



it I 



since one can see, through an identity similar to (|78D, that it has one doubly degen- 
erate eigenvalue. Then the integration measure on the matrix Q may be replaced by 
an integration over the eigenvalues ui,U2,t and a group integration, up to a Jacobian 
which may be found as follows. One considers a group element near the identity 

g = l + ie (84) 

with 

/ in n \ 

(85) 



ia b 
b 



with 



6) 



V 
V* 

Starting from a diagonal matrix Q it is transformed under this operation into 
Q + 6Q, with 



SQ = i[e,Q] 





v*(ui -U2) 



v{u2-ui) {it-ui)bii {it-ui)bi2\ 




-{it-ui)bii -{it-U2)b2i 
\-{it-Ui)bi2 -{it-U2)b22 



{it - U2)b2l 






{it - ^2)622 





(87) 
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The Jacobian J is then simply 

{it — UiY {it — U2Y ^ ^ 

The representation (^) leads to consider the supergroup equivalent of the HIZ 
integral defined by 

1 = j (i(^e^s*^^'5^"^ (89) 

where A and Q are diagonal matrices, A = diag(/i, /i, A, A) and Q = diag(Mi, U2, it, it). 
(We are still dealing here with a single ratio (^2|) of characteristic polynomials and we 
shall return to this point for higher correlation functions). The HIZ integral I may 
be obtained as the solution of a heat kernel differential equation. Indeed it satisfies 
the equation 

Aq/(Q) = eI{Q) (90) 

in which A is the Laplacian with respect to the matrix elements of Q, sum of second 
derivatives with respect to the five commuting variables, and to the four Grassman- 
nian ones, and 

Since the integral / is invariant under supergroup transformations of Q, it is a function 
of the three distinct eigenvalues of Q Ui, U2, it and one may express the Laplacian as 
a second order differential operator with respect to those eigenvalues : 

(The factor 1/2 comes from the degeneracy of the eigenvalue it). Let h denote the 
square root of the Jacobian 

h = (93) 

and substitute in the differential equation I{Q) = x/h. Then, with the help of the 
identity 

2hdt^ dth' hf^^duo,^ du^h' 
17 



1^ : v^-^— -^v^ m) 

2 dt^ ^1 9m,2 2h dt^ h ^1 ^ ^ 

the first derivative terms of x cancel in tlie differential equation and one obtains, up 

to contact terms which have a vanishing contribution to the final result, 

Id'x . ^ d\ ' 1 1 1 ^_ 

2 dt^ + duj + ^ (.t - u^y ^2{u,- u,f^ ^'"^ 

If one substitutes 

X = e-2^^*^+^"^'^+^"2^Y/|ni-U2| (7 (96) 

(the square root factor comes from the degeneracies in the matrix A), the pole 1 / {ui — 
U2) in (^) is cancelled. Then the HIZ integral I becomes 

in which g satisfies the Laplacian (heat kernel) equation, 

c^-^T^dg r.^T dg dg Id'^g d'^g d'^g I , dg dg 



dt dui du2 2 dt^ dui^ du2^ ui — U2 dui du2 

Remarkably enough, the solution of this equation is simply, 

_ 1 1 

^ ~ ^ ~ 2N{it-u,){\- fi) ~ 2N{it-U2)iX- fi)' ^^^^ 

Collecting all factors from (P?]), together with the Jacobian ( |SS|) F/v(A,/i) may then 
be written as 

J {uiU2)~ i'^'t - Ui){tt - U2) 2iV(A-/i) tt-Ui lt-U2 

^ ^-2imX+N(m+u,)^.^Nt--fiul+ul)^^^^^^^^ ^^QQ) 

From there one derives the density of states 
p(A) = ^Imlhn|-F^(A,/i) 

— -Im lim f dtdu- ^' ^ "^2! ^ , 1 



Stt'^N ij.^xJ (uiU2)~ {"^t — ui){it — U2) it — ui it — U2 

^ ^-2^Nt\+N{u^+U2)^l-Nt^-^^(ul+ul) {1^1) 
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This result has been also obtained recently by Guhr and Kohler[|l3]. Since this 



expression is quite different from (|68|) , derived in the section two, it is worth showing 
that it does gives the same result for the density of states. For simplicity let us 
consider the N=2 case. 



We return to the variables defined in (^) ; then p(A) reads 

p(A) = / }'^^' ^'^7^ ^e-^(*+^^)'-^(^~^)'~^''Mdp (102) 

' SttN J - p {{ity - 2itb + - pY y v i 

The imaginary part is evaluated as in (p§), 



^ SvriV J 'Sn^^^^^^ ^ ((^t)2 _ 2ttb + 62 - p)2^ ^^^''^^ 
1 ^ sgn(b)-!^^e-^(*+^^)'-^^'-^(^-^)'M (103) 



87riVy ' '(it -26)2 

The factor it — b in the integrand, is equivalent to taking a derivative with respect 
to A of the exponent. We also introduce an auxiliary integral over a variable s to 
express the denominator, 

2dX J-oo Jo Jo 

/CxD /'O /"OO 

-oo J —oo Jo 

This difference of integrals is an even function of A, since it may be expressed as, 

P=^|^[/(A)-/(-A)] (105) 

with /, after integration over t, is 

/(A) = r rrfsrf6se-^^(''-t+^)'+ 
Jo Jo 



2 I JV , 2 



dsse^-^^' I dbe-'^^^ (106) 







oo 



2 2JV 



Therefore, one obtains 
d 



aA JO 

= 2Ne-^^' + NXVaN^ +2N\ [ e'^dz (107) 

J-2NX 
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Dropping the odd terms in A, we recover the expression (^) for the N = 2 density 
of states. For general N, the same method apphes. Taking the imaginary part of 
(j)2_p)(iv/2) ; gives derivatives of the delta- function, and using integration by parts, we 
obtain the expression for the density of states 

5 Large N limit : the semi-circle law 

We now consider the large N limit of the density of states, based on the representation 



(|101| ). In the large N limit, one may use the saddle point method for the integration 



over t and Ui,U2. Since the denominators in ( |1(J1| ) have double poles, we again change 

itegral becomes 

it)^ 2it - 2b 



variables to 6 = "^1""^ , r = (l^—^faL^ Then the integral becomes 



2 ' 4 



/oo roo roo 

dt db dr 
-oo J~oo Jo 



M-A 7_oo J~oo Jo (62 - r)^/2 [(it)2 - 26(it) + 62 _ r]2 
^ ^-N{t+i\f-NbH2Nb,i-Nr (^]^gg^) 

Integrating by parts over r, one is led to 

p(A) = hm Im r dt r 1 e-^(*+^^)'-^(^-'')' 

Svr^ AJ^A j_oo j-oo b it — b 

+ - — - lim Im / / / dtduidu2- 



'(UlM2)^/2 it — u-^ 2UiU2' 



(109) 

In the second term, we have used again the variables Ui and U2 instead of b and p. 

The saddle points are 

t± = ^ (110) 

/i ± v//i2 - 2 

b±=u± = (111) 

The two leading saddle-points are, (t+,6+) and (t^,b^) . Other choices, such as 
(t+, 6_) give an oscillatory behavior, which can be neglected in the large N limit. The 
first term in ( |109| ) is same as GUE. For the second term, one needs to evaluate the 
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Gaussian fluctuation around the saddle point since there is a factor \ui — U2\. One 
retains also the pair of saddle points {t^, ^2+) and ^2-). The fluctuation 



around the saddle points, give a factor 1 / \J f"{tc) = 1/ ~ (2-^1/4 ^"^^ 

integration over t. Noting that 

lU -b+ = it+ -u+ = {\- ^)(-^=k=) (112) 



we find that the factors l/yf"{tc) cancel, and the second term becomes same form 
as the first term. We obtain by adding by the two saddle points, in the limit A ~ /x. 



= ^ 1„„ = i^2^ (113) 

TtN M^A a — /i TT 

which is indeed the semi-circle law in the large N limit for the density of states. 

6 Ratio of two characteristic polynomials and two 
point correlation function: preliminaries 

We now consider the ratio of two characteristic polynomials, 

^(A.,A......,^< ^::;^;:^;^::;^-^j > (114) 

for the GOE ensemble. We introduce again Grassmann variables 61,62 and complex 
numbers Zi, Z2 to represent these determinants. 

(Al, A2, /il, /i2) 
N 2 

IT TT aadZ(^nd6Qi(id6QiQ^ 
a=l a=l 

y ^ g*JV[eia(Al<5a(,-Xa(,)eib+e2a(A2<5af)-^af))t'26+2j'a(/^l'5ab-^a6)2;i!)+2;5£,(M2<5af,--'i^a!))226] ^1X5) 

We integrate out the random matrix X by (^5|). Noting that iiY"^ and tryy'^ are 



now 



trF^ = -{djp){6p6^) + {zlzp){zlz^) + 2{e^zp){zp^) (116) 

tryy^ = {ejp){epe^) + (44)(z^^,) + 2{e^zi){zpe^), (iir) 
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where a,P = 1,2, and (OaO/s) = Y.a=i^aadi3a- As for the one-point correlation func- 
tion, one introduces auxihary fields, bi,b2,b[,b'2, and use the Gaussian integral 
representation of (|^). Then, one integrates out 9,9,z*,z. The result is then cast 
into the form. 



F^(Ai,A2,/ii,/i2) = /[Sdetg]-^e"|strQ^+strQA 
where Q is an 8 x 8 super-matrix. 



Q 



( 


bi 


-Vi 


-V3 


-V4 




iV2 










bi 


-K 




im 


irji 




iVs 








b2 






iV6 




iVs 








-vl 


b2 


im 


ir]5 




ir]7 




-ir]i 




-iV5 




-tb[ 





-Wi 


-W2 




im 


if]i 


iVe 


im 





-ib[ 


-W2 






-m 








wl 


W2 


-lb'2 















W*2 


Wi 





-tbj 



:ii8) 



(119) 



Again one needs to diagonalize Q by a super-group transformation, and find the 
Jacobian for the representation in terms of eigenvalues. This yields a factor (ti — 12)^ 
for the (b', w) block matrices of Q, as discussed in a previous article[0. The Jacobian 
for the Grassmannian part is again obtained by linearizing the group transformation 
near identity as in (|S7|). Then, in terms of the eigenvalues of Q, we have 



Fn{Xi, A2, /ii, 



{i^h)itt2)] 



N 



n 



4 

i<j 



Ui 



Uj\{ti 



X / dtidt2duidu2du^dui 



3 "J 



(120) 



where / is the HIZ integral 



dge 



NStrgQg~^A 



:i2ii 



It is sufficient to take diagonal matrices Q = diag(ui, M2, U3, U4, iti, 'jt25 ^^2)? and 
A : A = (/ii, /Ui, /i2, Ai2, Ai, Ai, A2, A2). Note that the eigenvalues of A are doubly 
degenerate. This HIZ integral is again computed by a taking the Laplacian with 
respect to matrix elements, a generalization of what appeared in our previous work 
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The Jacobian J in this case is the generahzation of ( ]88| ) , which is derived by (|8 
The Laplacian equation is 



j^Ukh^_^Shzl^, (122) 



The factor 1/2 comes from the degeneracy of t. We set again h = J^/^ and substitute 
ip = x/h. Then, the same identity ehminates the first derivatives of x- Since 



the eigenvalues of A are degenerate, we factor out y {ui — U2){us — M4), and we have 
the following HIZ integral, 

(^1 - Ili<j - Mj(Ai - A2)2(/Ui - /ia)^ 



X ^'(ui - U2)(m3 - ^4)^e-2*^*^^^-2*Af*^^^+^("^+"^)'^^+^("^+"^)'^^ + (perm.il24) 

where (perm.) means sum over permutations of the Aj's and /ij's. The differential 
equation becomes 



(9ti dt2 dui du2 du^ du4 

2 ^ dti'^ dua^ ui - U2 dui du2 - U4 dus du^ 

2 V 1 1 1 1 



(^1 - ^2)^ (iti - 2(ui - ^3)2 2{Ui - U4)2 2(^2 - Us) 



2 



The solution to this equation may be obtained by expanding in 1/N ; it contains 
several factors 

ri ^ nri 1/1 1 1 1 \ i 
9 = [1 - T^][l + TT7( + + + ) + ■■■] 

I\ti2 4A' Mis Ui4 U23 U24 

i=l a=l '^a 

1 1 1 , , 1 1 ^ , 1 1^,1 1 ^ 

+ TT7^((— + —)(— + —) + (— + — )(— + —) 

AN^ Til Ti2 ri3 ri4 r2i r22 r23 r24 
23 



,1 1 1 1,1 1 1 1, 
+ (— + — + — + —)(— + — + — + —)) 

Til Ti2 ri3 ri4 T21 T22 T23 T24, 
1,1 1,1 1,1 1 1 1 , 

- + —)(— + —)(— + — + — + —) 

SM-^ Til T12 ri3 Ti4 T21 T22 T23 T24 

,1 1,1 1,1 1 1 1, 
+ (— + —)(— + —)(— + — + — + —)) 

T21 T22 T23 T24, Til T12 ri3 ri4 

1/1 1 ^/ 1 1 1 1 X/ 1 1m /X 

+ T7^(— + —)(— + —)(— + —)(— + —)] (126) 
IbA/^ Til ri2 Tis ri4 r2i r22 r23 r24 

where we define 

tl2 = {h - t2){Xl - X2) 

Ul3 = (Ui - U3){fli - ^12), Ui4 = {Ui - U4){fii - fi2) 

U23 = {U2 - Us){j2i - jj,2), U2A = {U2 - Ua){^Ii - ^12) 

Til = {iti-Ui){\i- Hi), ri3 = (iti - U3)(Ai - /i2) 

ri2 = {iti - U2){\i - fil), etc. (127) 

Note that all these variables Nt, Nt, Nu are of order one in Dyson's scaling limit. 
The series in 1/r terminate at order However the second factor in (|126|) is an 



infinite series in powers of 1/A^, and the scaling limit requires a determination of the 
full series. Remarkably enough a closed solution may be found in the scaling limit. 
This series turns out to be simply the HIZ integral for the GOE ensemble, and it will 
be discussed in the next section before returning to the two-point function. 

7 Heat kernel equation for non-unitary Itzykson- 
Zuber integrals 

Let us consider the integral, 

,NT:r(gAg-^B) Q28) 



in which g runs over the orthogonal group 0{2k), with the usual Haar measure of 
integration. We do consider group integration in this section, not supergroups . 
We shall also later generalize it to other values of /? with, as usual, P = 1,2,4 for 
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the GOE,GUE and GSE ensembles respectively. The eigenvalues of A are denoted 
• ■ ■ , U2k), and the eigenvalues of B as (/ii, ■ ■ ■ , fi2k)- Considered as a function of 
the matrix A the integral I satisfies a Laplacian equation, 

AaI = el (129) 

where is the Laplacian with respect to the matrix elements of A, and 

e = N^TTB^ = N^f^fil. (130) 



Since / is a function of the eigenvalues Ui of A one may replaced the Laplacian on 
the matrix elements by a differential operator on the eigenvalues 

1 (9 d 

where the measure J is given by the absolute value of the Vandermonde determinant, 

2k 

J=Y[\U^-Uj\ (132) 

i>j 

Let us investigate the simplest case of a 2 x 2 matrix {k = 1), for which one obtains 

(9^9^ Id Id 

[tT^ + TT^ + ^^TT ^TT^]^ = ^^33) 

awi^ 0U2 ui — U2 oui ui — U2 0U2 

To eliminate the first order derivatives, one substitutes 

^ X (134) 



V^l - U2 

Then ( |129| ) becomes 



""'^ + 1^ + i^-^X = (135) 



dui^ du2^ 2 {ui — U2Y 
One then substitutes further 

Noting that e = A^^(/ii^ + /i2^), we obtain 

df df (9V 5V 1 
2Nf,,^ + 2iV/i2^ + 1^ + 1^ + ^/ = 0- (137) 

OUi OU2 OMl C'M2 2,[Ui — U2) 
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From this differential equation one generates the expansion of the function / in powers 
of 

1 9 1 

^ ^ ^ + ^(^^I^^^I^Ki^;^ + ?,2N\u^-U2n^^l-^l2f ^ ^^N^^ ^^^^^ 

This series is in fact the expansion of a modified Bessel function, as may be 
recognized on the differential equation. Indeed if we look for a solution /, which is a 
function of the single scaling variable 

x = N{u^-U2){^ll- ^l2) (139) 
one sees that f{x) should be the solution of the differential equation 



fix) + fix) + ^/ = 0. (140) 



Therefore the k = 1 solution is 



/ = e^("^+"^)('^^+^^)/o(y (wi - W2)(/ii - /i2)) (141) 

where Io{z) is a modified Bessel function, whose asymptotic expansion is 

1 9 225 

in agreement with ( p.38| ) . 

In the expansion ( |138| ), we have implicitely assumed that {ui — M2) is of order one 
in the Dyson limit, where N{^i — ^2) is of order one. However if {ui —U2) ~ 0{1/N), 
the previous expansion is not valid and one has to perform a different expansion of 
the solution of the differential equation ( |137| ). Instead of ( |138| ) one expands 

/ = C^\u^ - n2|[l - - /i2)(Mi - M2) + ■ ■ ■] (143) 

in order to be able to deal with the regime in which the difference of the two arguments 
Ui and U2 becomes order of 0{1/N). It corresponds of course to the small x-expansion 
of the modified Bessel function (|142|) . 
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Let us now proceed to the k=2 case, which is relevant to our problem. We have 
J — ni<i<j<4 Wi — and consider the region ui > U2> > u^. The differential 



equation for I reads 



4 02 4 Id 

1:^.1 +Y.i:-^.^'=^' (144) 



1 

(i=l,...4), one has 



Substituting I — , ^ y, where A(u) is the Vandermonde determinant of the m 



A 92 1 „ 1 , , 

Again one substitutes x = e^"i''i+^"2M2+iv«3M3+iv«4/^4y^ and finds 

E + E #4 + ^ E ^^^/ - 0- (146) 

^ 97/^2 2 ^ (Mi - Mj)2 

This equation has beautiful properties in the scaling limit of interest. Before dis- 
cussing this propcrty,let us write the generalization of this equation for arbitrary (3 
{(5 = 1,2,4 are the standard values for GOE, GUE and GSE respectively), and for 
k- variables instead of four : 

^ Bf ^ B'^f 1 



in which y stands for 



?/ = /3(f-l)- (148) 



It is a remarkable property that this equation has a solution which is a function of 
the scaling variables 

Tij = N{iii - iij) {ui - Uj) (149) 

alone, and not separately of the u^s and the /x's. Given the origin of this equation, 
namely the HIZ integrals, it is clear that the solution involves only dimensionless 
products of the type /ij-Wj, and that it is unchanged by a simultaneous permutation of 
Ui and Uj accompanied by the same permutation on /li <-> /ij. However a direct proof 
or verification on the equation itself, turns out to lead to very elaborate combinatorial 
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identities if one tries, for instance, to expand the solution for large r's. The lowest 
orders are simple , but the calculations become very tedious, and non trivial, at 
higher orders. Let us reproduce here simply the first terms : 

i<j ^'tj i<j 'ij ° i<j '*i 

This expansion shows our point : the coefficients are pure numbers, completely in- 
dependent of the parameters /i's, once the function is expressed in terms of the Xj/s. 

This is a generalization of the expansion of the modified Bessel function( [L38| ) ; 
one recovers the simple Bessel limit if one lets all the m's, except the first two,go to 
infinity. 

As discussed hereabove, all we need in the large N limit, is a solution for the 
degenerate case, /^i = /i2 = and /i3 = /i4 = /i'. To deal with this degenerate case, 
one writes / as 

/ = \J\{Ul- U2){U3 - U4)\g (151) 

From (|146|) , we have 



OAT ^9 , o Ar ^9 . ^.r ,dg ,dg d^g 



dui du2 dus dui ^ dua 

+ ^ I ^9 _ dg_ ^ I ,dg dg 



2 



Ml — U2 dui du2 Us — U4 dus dUi 

,11 11 11 11, , , 

2 [ui - uzY 2 [ui - u^y 2 [u2 - u-iY 2 [u2 - u^Y 

This function g may then be obtained as an expansion in powers of Terms 
involving poles in — ^ — and — ^ — do not appear in the expression for q. We have 

^ 1 1 1 

1 

^ 4Ar(M2-M4)(/i-yU') 

9,1 1 1 1 , 

+ :::rTTT77 ^7 ttt + 7 ttt + 7 ttt + 



32A^2(/i - (mi - M3)2 (ui-U^Y {^2 - U^Y {^2 - U4) 



2 J 



1 r 1 1 
+ ^7 r + 
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3 1 1 1 

16N^{fX - IJ,')'^ (Ui - Us){ui - Ui) {Ui - U3){U2 - Us) {Ui - U4){U2 - U4) 

+ 1 \-i (153) 

This expression has a comphcated structure ; one may write the various terms 
as diagrams, but the coefficients depend upon the topological character of those 
diagrams. For instance, the coefficients at order 1/A^^ the three types of diagrams 
have weights 9/32, 1/16 and 3/16. However, we find that in the large N limit, the 
expressions simplify. 

In the large N limit, if we concentrate on the saddle points u\ = = and 
M3 = = M_, one obtains 

1 2 1 

^ " ^ N{u^ - u^){n - n>) ^ m{u+ - u^yifi - n'f ^ *^^7V3^- ^^^^^ 

If we introduce the scaling variable 

X ^ N(u+-u^)(Ijl- 11'), (155) 
one finds that g{x) satisfies the differential equation. 

The expansion of g in powers of l/x is 



9=E5^ (157) 



n=0 ^ 



It is interesting to generalize this formula to general k, {k > 4). Denoting Vj — 
tifc+j,(j = l,...k), one deals with the differential equation. 



2 



+ 1 ) \ Y ^ ( ^^ ) 

+ ^EE[^7-^]=0 (158) 

29 



For general k, the situation is similar. In the scaling limit, in which the leading 
saddle points involve only two distinct values of the form Ua = m+ and Va = u^, one 
introduces the same scaling variable x (|155|) , and obtain the differential equation, 

dx'^ ~'~ X ^ dx ~^ Ax"^^ ^ ^ 

An expansion in powers of 1/a; of the solution follows easily : 

g|M». + 2)(t + 4)...(t + 2(p-l)ff 

One recovers the previous results for A; = 1 and k = 2 from this expression. 

We have considered hereabove the GOE measure </ = JJ^ — Mj| • Let us now 

i>j 

go to arbitrary /5, i.e. work with J = JJ^ l^j — Uj]'^ ; the GUE integral corresponds 

i>j 

to (3 = 2 and the GSE to (3 = 4. When the eigenvalues /ij, which are set at the 
saddle-points in the scaling limit, degenerate again into two groups, one equal to /i 
and the other one to fi', one finds through identical steps 

OAT \- ^3 ,^ dg ' d'g ' d'g 

a=l ^'^a a=l Q a=l "-^""a a=l ct 

1 , dg dg . ^ ^ 1 . dg dg . 

^pUo,-Up dUa dUp o^f3'"a-Vp dv^ dVp 

In the scaling limit, in which the m's and the w 's approach repectively m+ and u_, one 
introduces again the scaling variable (|155|) and one obtains the differential equation 

9^9 k-1 dg A;2/3(2-/5) _^ 

dx'^ X dx ^ 

which generalizes the previous GOE expression for /3 = 1. The expansion of the 
general /?, arbitrary k solution, in powers of 1/x follows 

^ , ^ (/?fc)((2 - miPk + 2)(2 + k{2 -P))... {2{p - 1) + k{2 - 
^ = 1 + E (163) 
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In the unitary case, (3 = 2, the solution reduces to the first term g = 1 : this is of 
course the well-known Itzykson-Zuber result which is semi-classically exact. 

For /3 = 4, and = 1, we also see that the asymptotic expansion ( |163| ) stops at 
first order, since the terms with p > 1 vanish in the series ( |163| ), a fact that we had 
already found and used in ( 126 ) for integrating over the tj's. For /3 = 4, and = 2, it 
stops at third order. The fact that the Itzykson-Zuber integral for this scaling limit 
of the GSE, is semi-classical with only a finite number of corrections had already 
been discussed and used in 

Therefore one sees that, in the scaling limit of interest, in which we deal with an 
Itzykson-Zuber integral for a degenerate case, there is a remarkably simple expression 
for arbitrary (3 and arbitrary dimension k of the group integration. The expression 
is either semi-classically exact for the GUE, corrected by a finite number of terms 
for the GSE, or an infinite, but explicit, series for the GOE. The fact that k , the 
dimension of the integral, appears as a parameter will allow us later to continue in k 
and use the replica method in the k goes to zero limit. 

Returning to the GOE, let us note that if we had included an i in the exponent 



of the Itzykson-Zuber integral in ( |128| ), and multiplied by the factor e*^/x, we would 
have obtained for k = 2 

.e^^ , , cosx sinx 2cosx 6sinx 
Re — g[ix)\ = + — ^ — + ■■■ 

^ (fc-l)!sin(x + ffc)(-l)(^-^) 
k=i 

which is the large x asymptotic expansion of the integral 

r°° sin z 



J 

Jx 



dz. (165) 



This integral appears in the scaling limit of the resolvent-resolvent correlation func- 
tion fo the GOE ensemble. 
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8 Two point correlation function 



We are now in position to return to the two point correlation function of the GOE, 
which will be deduced from the ratio of products ot two characteristic polynomials. 
Indeed one may obtain the resolvent-resolvent correlation function 



1 d d det(Ai -X)det(A2 -X) 
MAi, \2) - ^2Qx,dX; < det(/.i-X)det(/.2-X) ^ ''^-^-'^-^^ ^^^^^ 

whose double discontinuity gives the two-level correlation function. For this correla- 
tion function, the relevant terms in the scries g discussed in the previous section, are 
those which involve a factor (— -l- — )(— -l- — ) since this factor yields pole terms 

^ Til ri2 ' ^ T23 T24 ' J r 

of the form tt h: r. Therefore the relevant part of g, after taking derivatives 

(Al-/il)(A2-M2) S:' HI O 

and focusing on terms with poles at /li = Xi, is 
1 

4/V2' 

1 1 _ 1 1 



9l - -rvTKgtgu 



'iti-Ui iti-U2 it2-U3 it2-U4 (Ai - //i)(A2 - //2) 

— )- 

2—U2 A2 

(167) 



x(l ^("^ + '^)^)(1 "^("^ + "^)"^) 

2N iti — Us iti — U4 Ai — 112 it2 — Ui it2 — U2 A2 — Hi 



with 



ft = 1 - ^ (168) 

In the large N limit, we have three types of saddle points, i) {t^ , t2, uf, U2, uf, uf) 

or {ti,t2,u^,U2,Us,u^) ii) {tf ,ut ,u^), {tt,t2,u];,U2,ut,ut), 

{t];,tt,ut,ui,u^,u^), or {t];,tt,u];,U2,ut,ut). iii) {tf ,t2 ,ut ,U2 ,ut ,u^), and 
similar combinations. 

In those expressions the saddle points are given 

t- = Z!A±^ (169) 

= ^^^^f^ (170) 
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with Ai, A2,/ii,/i2- 

For the saddle-points of type i), is simply 



9u = \J\{ui- U^){ui - Ui){u2 - Us){u2 - U4)\{j2i - 122^ (171) 

Thus we have, as partial contribution of this saddle-point to the correlation function, 
p»(Ai,A2) from (|124D and (|T67|) , 

(Ai-;Ul)(Ai-/i2)^(A2-yUi)2(A2-Ai2)(— ^ + — ^ )( ^ ^ 



iti — ui iti — U2 it2 — U3 it2 — U4^' 
r i 11 11 

~ N{ti-t2)iXi- X2) ~ 2NHti -U3^ ih -uj\^-^i2 

[1 ^ ( ^ + ^ ) ^ ] ^ 

2N it2 - ui it2 -U2 A2 - Ail (Ai - A2)^ 



Af(tf+t|)"^(u2+„2_,_„2_,_„2)^^(„^_l_„2)^j+Ar(M3+u4)^2~2iAftiAi-2jAfi2A2 



+ perm.(Ai ^ A2,/ii ^ /i2) (172) 
from which follows 

W^Ai,M2^A2C'A2C'Ai J {uiU2U3Ua)~ Ilj=l Ila=l{'^'tj - ^a) 

(Ai - ;Ul)(Ai - ;U2)(A2 - /Ui)(A2 - /i2)(— ^ )(t 



iti — Ui iti — U2 it2 — U3 it2 — Ui' 
1.1 1 1 1.1 

+ )(^ 



AN"^ iti — Us iti — U4 it2 — ui it2 — U2 (Ai — A 



2 



\2 



^-N(tf+t^)-f(uf+u^+u^+ul)+N{ui+U2)fii+N{u3+U4)fi2-2iNtiXi-2iNt2X2 (173) 

For those saddle points of type i), the difference ti — t2 is of order 1/A^. Then, 
cancelling terms between the numerator and the denominator, one obtains 

p^^Ai,A2) ~ \m dtdu- -jr— rr- r(- + - ) 

J {uiU2U3U4)~ [^tl - Ui)[lti - U2) tti - Ui iti - U2 

1 , I 1.1 



(ih - u.3){it2 - U4) it2-us it2-U4 N"^ 

g-Af(ti+t2)-2jA'tiAi-2iArt2A2-^(«f+«i+«§+«|)+Af(Mi+M2)Ai+Af(u3+M4)A2^]^y^^ 
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which is just a product of density of states as given in the equation ( |101| ). Therefore, 

by considering the saddle point of type i) (tf, uf, U2, u^, M4 ) and (t]^ , u{ , U2 , , u^) , 

we obtain 

p»(Ai,A2)=p(Ai)p(A2) (175) 

We now return to the integral ( |173| ), and consider a saddle point of type ii) ; one 
thereby obtains the sin- kernel squared, since now terms like {iti — ui) are replaced 
by the density of states p(A). The exchange of pairs between (ui,M2) and (^3, U4) in 
the denominator gives 



This term adds to (|175| ), and by adding both terms, one obtains the expected van- 



ishing contribution in the limit Ai A2. 

Now, we take into account the remaining terms, gt and Qu- In the case ii) ui — U3 
is now proportional to the density of states p, which is of order one. Then all terms in 
Qu are of the same order in the Dyson limit, and we have to return to the discussion of 
the previous section. For the saddle point of type iii), such as (tf, t2 , uf , U2 , , u^) , 
one has 



p(Ai,A2) = Im J dtdu 



(iti)^(it2)^ 2iti — Us — U4 2it2 — Us — Ui 



)- (zti - usY{iti - u^Y {it2 - uiY{it2 - U2) 



Nnp{\^ - A2) ^ ^ 4iV7rp(Ai - A2) ^ " 

^ ^-N{tl+tl)-2iNti\i-2iNt2\2-f(ul+u^+ul+ul)+N{ui+U2)Xi+N{u3+U4,)X2 (^177) 

in which we have replaced various terms by p(A) and (Ai — A2) by using their saddle 



point values. We note that tf — t2 = \/2 — }? = vrp(A). One introduces the same 



integration variables b and p already used in (|67|) , and integrate over pi and p2 as in 
( 109 ). Then, there is a part of the integrand which has a double pole in ^^^_^^^^2 in 



the large N limit. The imaginary parts are then taken for Ai and A2 independently. 
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Thus one obtains for gt , after multiplication by a factor 



iVp(A)(Ai-A2)' 

N7rp{\i-\2)^ iV7rp(Ai - Aa)^ 

cos(7rA^p(Ai - A2)) sin(7rA^p(Ai - A2)) 

Arp(Ai - A2) NVp^Xi - A2)2 
d / sin X ' 
da; V X 



(178) 



where x stands for x = A^7rp(Ai — A2). Next we return to the Itzykson-Zuber factor 

Qu of (11531). After multiplication of Qu by -— , one obtains 

A/p(A)(Ai - A2) 

cosx sinx cosx sinx 
= + ^-2^-6^ + --- 

lAj iXj iXj 

dz (179) 

z 

Therefore, we have by multiplying these two factors, 

, , d ,sinx, r°° smz , , , 

The third type of saddle points iii), for instance {ti ^ul ,U2 ,u^) , does not 
yield any imaginary part, and may thus be dropped in this GOE calculation. However 
for GSE, they do contribute as well. 

Adding the type i) ( |175|) and ii) ( |176| ) results, one obtains the two-point correlation 



function for the GOE in the large N limit. 



9 / , N r , sin X , , (i , sin X , /'°° sin 2: , ^ 
p(Ai,A2)=/(A)[l-(^)^--(^)xy^ ^dz] (181) 

where x = 7rA^p(A)(Ai — A2). Comment : This result is of course well-known ; it 
has been obtained long ago through the technique of skew orthogonal polynomials. 
The point of the long derivation presented here, through generalized Itzykson-Zuber 
integrals, is that it can be repeated in other cases, such as non-invariant measures 
involving an external source, for which the standard method does not apply. 
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9 Zero replica limit 



Up to now we have computed averages of ratios of characteristic polynomials defined 
in (|l|). The usual two-point correlation function of the resolvent operator is obtained 
from this ratio by differentiation, as shown in ( |166D . Instead of such ratios one may 
use an alternative "replica" method, as follows. 
Consider the correlation function 

F2k{Xi, A2) =< [det(Ai - X)]*^[det(A2 - X)f > . (182) 

Since [det(A — X)]'' = exp[A;tr log(A — X)], one may recover the correlation functions 
of the resolvent, by letting the replica number k go to zero : k —>■ since : 

Let us first discuss the zero-replica limit ( |182| ) for the GUE, and show that the 
well-known exact expression for the two-point correlation function of the resolvent, 
is correctly recovered. There were discussions of the validity of the replica limit in 

In an earlier work lQ,the universality of , in the Dyson short distance limit 
was proven and its expression was found to be (for a probability measure proportional 
to exp-NTrV{X)) , 

F2,(Ai,A2) = e^(^(^^)+^(^^))e-^^i(2vrp(A))^^ 

kl 

n ^exv - MY u ) ^'^^^'-'^^^ (184) 

a special case of the general formula derived in for different Aq,'s : 

2k 

i^2fe(Ai, ■ ■ ■ , A2fc) =< n det(A« - X) > 
0=1 

(185) 
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in which one has used the scahng variable x = A^7rp(A)(Ai — A2), xi = 2ttNp{X){Xi — X) 
and A = (Ai + A2)/2. This contour integral is a compact representation of the sum 
over the (2A;)!/(A;!)^ saddle points which govern the Dyson limit, and it is particularly 
useful for degenerate cases in which some of the Aj's are equal. Indeed in those cases, 
subtle corrections are present in the sum over saddle points, which are not easy to 
handle if one lets the A's approach each other too soon. The poles due the degeneracy 
of the Aj's have to be cancelled by the sum over permutations over all possible saddle 
points. The formulae ( |184| ) and (|185|) do give a correct answer for the degenerate 
cases of equal x/'s. 

The calculation is done as follows. One first shifts — Mq, = f „ + x and drop 
all factors in F2k, which approach one in the zero-replica limit. Then we have to deal 
with the integral 

^ = / n ^ -p(-*- - ' E n.„Si +'^/«.)' 

which vanishes if one sets A; = in the integrand. The zero-fc limit is slightly tricky 
since the number of integrations is also equal to k, and one has to find explicit 
answers for the integrals before approaching the required limit. We are interested 
here in terms of order k"^ for small k. Thus one expands the denominator in powers 
of X, 

^r^-l-4^ + 0(.^) (187) 

Since the one gives a vanishing contour integral, up to order /c^ it is sufficient to 
expand one of the factors (1 + 2x/v)^^ to order k, since there are k ways of singling 
out a particular Va ■ If one expanded both (1 + x/f^)"^^ and (1 + x/v/^)'"^^, the 
number of choices k{k — l)/2, together with the two powers of k, would give a term 
proportional to k^. 

Then , in the zero replica limit, it is sufficient to examine the integral 

/(I) (a;) = ke-^'^ /rr — e-^^^'-T "^'^^j , (188) 

^ ' J 271 (l + 2x/t;i)^n?t^P 
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This integral makes it clear that there are oscillatory terms of the form e^^^, the 
contribution of the pole Vi = —2x, and a non-oscillating term which comes from the 

Sill cc 

pole Vi = 0. We will show indeed there are such terms and they are ( One 

X 

then expands e~^''^{vi + 2x)~'' in I^^\x) in powers of x. This generates the following 
integrals, 

( (2+P)(4 + p)---((2fc-2)+p)(p-3)(p-5)---(p-(2fc-l)) 

/=o(^ + 0! (2A;-l+p)! 

^ z=o(^ + 0!r(2A: + p)r(l + f)r(A:-f)r(3-p) ^ ^ 

Note that when p is an odd integer, 7^^'' does not contribute to /^^-'(x), since one has 
to take a real part. When p = 0, "yj?^ coincides with a well-known universal number, 
related to the moments of the Riemann zeta-f unction. In the zero replica limit of 
A; — > 0, we have found in an earlier work that that 

l'i3n(^ = l + *=(H--)H-0(*=) (190) 

where c is Euler's constant. Therefore we have 



By expanding (^^^^2x)'' powers of x, we obtain 



00 7,2 1 

/«(a:) = E Hn-){2xr- --— - (192) 

p=even P {^-pW{P+^) 

Thus the second derivative of J*-^-* with respect to x is finally given by 

q2 r(l) 00 r)p 2ix I -2ia; _ 9 

^ = -PRe S i-fAx"-' = -^'^^I^S ^ 

L/o. n=even F- 



2^2^sinx^2_ (193) 
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Finally there is an additional constant term which comes from the second derivative 
of e"*'^^ in J, which was neglected in Therefore adding this constant, we obtain 
the well known two-point correlation function of the GUE by this replica method, 

p{X„X2)=p\X)[l-C^r] (194) 

X 

We may now we proceed to the GOE case. The Itzykson-Zuber integral for the 
GOE case has been discussed in the previous section, when we dealt with ratio of 
characteristic polynomials. We now consider the following moment 

^ [det(Ai - X)]'=[det(A2 - X)]^' ^ ^^^^^ 
Again the zero-replica limit will be used to obtain the two-point correlation function 
of the resolvent operator for the GOE. 

Indeed again here we deal again with the Itzykson-Zuber integral for a degenerate 
source, with only two distinct eigenvalues, Ai and A2, both fc-times degenerate. In 
such cases the heat-kernel satisfies the differential equation (|159|) . The asymptotic 
solution of the solution in the scaling limit, is given by (|160| ), 

I f. [k{k + 2){k + 4)...{k + 2{p-l))]' 

In the zero-replica limit k —>■ 0, we need to keep only the terms proportional to k'^ . 
To order fc^, we have from ( |196[ ), 

A;2 (p-l)! 
4 p pxP 

Then taking the second derivative of gfwith respect to x, we have a factor {p\ + {p — 
l)!)/x^^^, which is precisely what one obtains from the asymptotic expansion of the 
standard two-point correlation function of the GOE ensemble , 

p,{x) = p2(A)[i - - r 'J^dz] (198) 

X (XX X Jx z 

Indeed taking two derivatives with respect to x of the second and third terms : 

\ dp , r°° sinz , .sinx,, d ,sinx, r°° smz , 

-— / dzf = i ?-^{ / dz 199 

2dx^^Jx z ' ^ X ' dx^ X ' Jx z ^ ' 
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Next let us compare our expression of g with tlie factor ( / dz^. /^From the 

Jx Z 

asymptotic expansion of the integral 

(200) 

one obtains for the square of this quantity 



1 -2ix.V- P- 



— 



'"^T.7£-,m' (201) 



4 2xP+i 



Indeed our expansion ( [197| ) for g in the zero-replica limit agrees with the first term, 
the non-oscillating term of the above equation, since 



The Itzykson-Zuber integral which led to ( |196|) , has been derived dby a saddle- 
point method , choosing uj = m+, (j = 1, k) and ui = u^,{l = k; 1, 2k). This 
gave a factor e"^^^ which reduces to unity in the limit A; ^ 0. Those saddle-points 
contribute to the non-oscillating terms. However there are other saddle-points that 
one needs to consider in order to recover the oscillating part. This is analogous to 
a phenomenon recently analyzed by Kamenev and Mezard As for the GUE, 

the zero-replica limit requires only the e^^*^ oscillating terms. Those terms may be 
obtained through saddle points of the following type : one divides the Ui (i = l,...2k) 
into two groups. M2, Mfc) and {uk+i,Uk+2, ■■■,U2k)- One chooses one Ui to be 
a u_ in the first group, and the others u+ ; similarly one u+ in the second group, 
and the others are m_. For instance, in the first group, ui = u^,U2 = u^^u^ = 
u+,...Uk = M+, and for the second group u^+i = u+,Uk+2 = u_,....,U2k = u_. The 
combinatorial factor summing over all such choices is k'^. The differential equation 



for this degenerate case is similar to ( |152|) , but the combinations of {u2 — Uj) for 
j = k + 2, 2k, and (wj — Uk+i) for i = 1,3, k are eliminated. One then modifies 



40 



(|T5lD to be 

/ = \J\{ui- U2){U3 - U4){Ui - U^){U2 - Ui)\g (203) 

for this purpose, in the case k = 2 given as an exmaple. Then g satifies a differential 
equation which yields the terms proportional to e^*^ as a powers series in 1/x up to 
order The first order is 

= fM^H^ (204) 

and it vanishes in the zero- replica limit. At second order, one has 

(2) _ 9A;2[(A;-1)2 + 1] ^ 3e{k-lf{k-2) , ^^[(fc - - 2)^ - 2(A; - l)^^ 



32x2 16x2 32x2 

(205) 

which come from the contribution of three different diagrams connecting two lines; 
a double line, a connected line, a separate two lines. The diagrams are the same as 
for the non-oscillating calculation, but the weight factors are different. In the zero- 
replica limit, it gives 5'*^2) = k"^ /Ax^. At third order, we have 7 different diagrams. 
Adding thir contributions, one finds in the zero-replica limit. 

If one compares this series with the terms proportional to e^^^ in ( |201|) , 

4^^ 2xP ^ ^ ^ Ux2 2x3+ J 

one sees that the two agree up to this order. 

We have a term ^-73-. By the second derivative of this term, we also have A;2 term 
x'^ 

as 

(i2 e^''^ , n , 1 



k\l--) + Oik') (208) 



(ix2 X*^^ x2 

^—ikx+2ix 

We also have similar terms from the derivative of — . Together, we have 

x'^ 

Sill X 

1 — 2( v. Adding this term to (|199D , we obtain the two-point correlation function 

X 

of GOE, 

9/,xr /sinxX^ d sinx f°° sinz , , , ^ 

MX,,) = /(A)|l - (— ) - ^— / —dz] (209) 
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Thus we have shown here that rephca hmit of the moment of the characteristic poly- 
nomial gives the consistent result with the well known resolvent two-point correlation 
functions both for GUE and GOE. 



10 Gaussian symplectic ensemble 

The Gaussian symplectic ensemble (GSE) is easily formulated as an extension of the 
GOE. Let X be a quaternion symmetric matrix. Let us consider as an example the 
N=2 case ; Xu and X22 are both real numbers , whereas the off-diagonal element 
X12 is a quaternion and X21 its conjugate. The quaternion X12 may be written as 

Xi2 = a + bei + ce2 + de^ (210) 

where Cj are the quaternion basis (i.e. up to a relabelling and a factor the Pauli 
matrices ) ; in the basis e, = C O.Ve. = '}] = ( the 



X' 



(211) 



,0 -zy ' " VI ; v-^ . 

coefficients a, b, c and d are real numbers. One can write instead the matrix X as 
ordinary 4x4 matrix X', the elements are then usual complex numbers, 

/ xii u iv \ 

xii iv* u* 

u* —iv X22 

\-iv* u X22/ 

where xu and X22 are real, u and v complex. The relation [det(A— X)]^ = det(A— X'), 
aloows us to write 

[det(/.-X)]2 det(/.-X') > ■ ^2^2) 

The density of state p(A) may be easily deduced from F2. From the relation ( pi2| ), 
one obtains in the N=2 case, 

^-'^■^'-// i^:::!'!;!:g <'--'->'--''-"''lrlr f^^^' 

in which ti and ^2 are the eigenvalues of X. The imaginary part of F2, from which 
one deduces the density of eigenvalues through 
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is obtained by picking up the contribution of the pole ti = jj,. Taking the imaginary 
part, a derivative with respect to A and setting /i ^ X afterwards, one obtains 

which may easily be checked directly. 

For the inverse of the characteristic polynomial, we find easily the following for- 
mula in the X' representation (still for the N — 2 example) , 



1 

det(/i-X') ^ 



•' i=i 
i=i 
^ •' i=i 

The imaginary part is the Hermite polynomial H^[ij). 

For general N, similarly, we find that the expectation value of the inverse of the 
characteristic polynomial is 



and its imaginary part is simply H2n-i{iJ')- 

For the inverse of the product of two characteristic polynomials, one writes 

(det(//i-X))2(det(//2-^))' i [det(5)]2^ ^ ' 

where X is an A'^ x A'^ quaternion symmetric matrix, B a 2 x 2 quaternion matrix, 
and M is a diagonal matrix M — diag{iii, //2)- 

The average of the square of the characteristic polynomials may then be written 

as 

< [det(A - >= J dA[det{X - A)fe-^''^' (219) 
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where A is a 2 x 2 real symmetrix matrix. The quantity [det(A — A)]^ is a polynomial 
in A of order A^^. 

Finally the ratio F/v(A,/i) may be written as an integral over a super-matrix Q , 
the derivation being similar to that for GOE, 

A super-group diagonalization and the Itzykson-Zuber integral ( for which we may 
use the same formulae as for GOE since the Jacobian has the same form after the 
diagonalization) leads then to 

^ X /■ , , , itit2)^ \ti - 12\(\ - a? , 1,1 1 

J u^^^ [u — lti)[u — lt2) N{X — fl) U — lti u — lt2 

(221) 

The density of state p(A) follows : 

p(A) = <lTr5(A-X)> 



lim Im / dtidt2du- 

J-oo 



(^1^2)^ \tl - t2\{2u - ih - lt2 



N 



-^{tl+tl)-Nu^-iN{ti+t2)X~2iNufi 



The formula is quite similar to that of the GOE, except that the combination ^ ] J'^ 
is raised here to the power 2N instead of — A^. This difference makes the calculation 
of the imaginary part easier, since the contour integral on b gives a contribution to 
the imaginary part, similar to that of GUE, and therefore the result does not involve 
the incomplete Gaussian integrals (such as B{x) in (|5P|)), which appear in the GOE 
case. 

Let us consider now a ratio of characteristic polynomials F/v(Ai, A2, /ii, /U2) defined 

as 

F (\ X n n) -< [det(Ai-X)det(A2-X)P 

F^(Ai, A2, /ii, /i2) -< [det(/.i-X)det(p2-X)]2 > ^^^^^ 

With the supermatrix formalism and its supergroup diagonalization, one can write 
again an integral over eigenvalues, similar to ( |120| ) in the GOE case, 

J (UiU2)2^ Y{i=iY{k=i{,Ua, - itk) 
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X Idtidt2dtsdt4,duidu2 (223) 
where / is the HIZ integral, 

^NSUgQg-^A (224) 



Jdge' 



withQ = diag{iti,it2,ih,itA,ui,ui,U2,U2), and A = dia5i(Ai, Ai, A2, A2, /xi, //i, //2, A*2)- 
The HIZ integral is in fact identical to that for the GOE case with corresponding 

variables. In the supergroup integration for the GOE case one found a product of 

two series, one infinite, and the other one finishing after a finite number of terms. 

The same structure appears in the present supergroup integration for the GSE case, 

except that the finite and infinite series have their resepctive variables switched. 
Thus in the large N limit, one has to consider, as in section 8, three types of 

saddle-points ; the saddle points of class i) and ii) give the same answer. The saddle 

points for ti and Ua become 

-a±V4^ //±iv/4^ 
t± = u± = (225) 

Note that for the present GSE case, the density of state p(A) is given by p(A) = 



V4 — A^/27r. Therefore, by the same arguments as for the GOE, we obtain the two- 
point correlation function 

/X > V /Sin2a;,r, d ,sin2x, f°° sm2z , , 

ft(A:, A.) =p^(A)[l-(—)^ --(—)/ —,4 (226) 

We have neglected the third class of the saddle point iii) for the GOE. In the 
GSE, those saddle points do contributite, since the imaginary part is taken from a 
contour integral over 61 and 62- Since the HIZ formula for /3 = 4, A; = 1 takes a simple 
form, we obtainfor the saddle point values of those b's {b^, 6^), and thus 

d sin 2a; 

/ = (227) 

dx 2x ^ ' 

For the corresponding t integral, we take the following saddle points : a set of saddle 
points such as {ti ,t2 itl) gives the non-oscillating constant contribution. Note 
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that ttt^ = 1, and iN{ti + )Ai = -A^A^. Thus there is a correction to ( |226| ) by 



(|227|) . By taking the normahzed coefficient, which makes the two-point correlation 
function becomes O(x^) in the small x limit , a property of the GSE, we find the 
scaling limit of the two-point correlation function, 

2,,.f, , sin 2^ 2 (i sin 2x , . /-"^ sin 22: tt \ 

= (1 - + /; (228) 

11 Extension to an external matrix source 

Let us consider an external matrix source A be coupled to a real symmetric random 
matrix, or to a quaternion self dual random matrix. The corresponding Gaussian 
probability measure 

p^(X) = 1 -f trX^+mrAX ^229) 

Zj 

has lost the invariance under GOE or GSE. 

In a previous article 0, we have considered this external source problem for the 
correlation functions (na=idet(Aa — -^)), in which X is a real symmetric random 
matrix. This is done by integrating over X 

where Y = — Y%t=\ ^aJ^ab- One may assume, without loss of generality, that A is a 
diagonal matrix. Then, the only new term with respect to the zero-source case, is 
exp[— iA^ ajOajOaj]- This diagonal term modifies the previous determinant (detS)^, 
and gives instead IljLi Il^=ii'ti ~ T'^-j)- In this way we had obtained in [|] that, when 
all the Aj are equal to A, one has 



k N k 

([det(A - X)r) = e-^i:^? / n IK^^ - Uiti - t.)V^i:*?+2^^^S*' n dti. 

■' 1=1 j=l Kl' 1=1 

(231) 
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Similarly for the ratio of two characteristic polynomials, we have Y = — Y,a=i (^aa^ab- 
Zaa^ab as in (^). The external source gives a diagonal shift for 9 and z. This leads 
to the modification of the Sdet term. For instance, in the k=l case, one replaces the 
super-determinant [Sdet]"^ by 

^ K^^^) (232) 

[SdetQ] 2 ^=1 [{a^ — Ui){a^ — M2)] 2 

Then the density of states, for the GOE modified by an external source matrix, is 
given by 



p{X) = -^ImhmJdtduU 



A'-Aj ^J^ [(a^ _ _ ^2)]! (it - ui){it - U2) 

it — Ui it — U2 

Let us consider the one point Green function G{z) (with an external source) , 

Giz) =< ti—^ > (234) 
z — X 

for a real symmetric random matrix X. It is given by the derivative of F(A,/i), 

^. . ,.9 det(A-X) 

G(z) = hm — < — ^ -f > 

^ ^ A=/.=^9A det(/i-X) 

a-y — it \ui — U2\{2it — Ui — U2) 



7=1 V""! - CL'y){u2 - a-y) - UiY{it - U2Y 

^ ^-2iNtX+N(u,+u,)t.-Nt^-l^[ul+ul)^^^^^^^^ (235) 

Using the change of variables, h = HJim^ ^ _ ("1-^2) ^ obtain 

/^f^^l _ f fjf f rlh n dr ~ -2iAftA+2iV6M-Aft'~Af(f'2+'-) 

^ n(52_^_2a,6 + a,2)i/2 (^36) 
After integration by parts over r, we obtain an expression, similar to the GUE case, 

G(-) = / n(i - 1™^) 

J ^ i\[b — ayj it 2m 2tt 
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In the large N limit, one may make the replacement 

Si, It ^ rt^ I , 

^=1 Nib -a^) N ^-^b- 

Let us denote the (non-random) density of states of the external matrix A as poia), 
we may then write the r.h.s. of (|238|) as 

exp(-zt / da^^^). (239) 
J u — a 

We define the resolvent of the external source Goyz) 

G,(z) = [ da^ (240) 
J z — a 



^From ( |237| ) and ( p38| ), we obtain 

dG r db 1 



(241) 



dz J 2m u + Gq{u) — z 
The contour surrounds all the eigenvalues a^. The zeros of the denominator satisfy 

M + — V = z 242 



As discussed in a previous paper ||15|, we take the poles of u{z) = z — ^ + C'(^) and 
^ = 00, and then 

^ = 1 ^ ^ Mzi 

dz 1 + ^^ dz ^ ' 

+ du{z) 

The integration gives 

G{z) = z- u{z) (244) 

Since u{z) is a solution of m + Go{u) = z, we obtain the following equation, due to 
Pastur |2§ for the GUE 

G{z) = Go{z - G{z)) (245) 

Thus we obtain the same Pastur equation as for the GUE case ; this is easily un- 
derstandable from a diagrammatic analysis ; in the large N limit, planar diagrams 
are simple rainbow diagrams and do not distinguish between GOE, GSE or GUE at 
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leading order in the large N limit. For the GSE, a similar algebra would lead to the 
same equation. 

For the two-point correlation function , k=2, the same shift for the Sdet occurs 
in presence of the external source. We have found in the previous sections that the 
resolvent two point correlation functions in GOE and GSE in the Dyson scaling limit. 
We now discuss the two point correlation function in the Dyson scaling limit when 
the external source matrix is coupled to the random matrix. We have already given 
a proof of the universality in GUE. The argument goes as follows for this GUE case 
PI- 

In the presence of a matrix source there is a kernel Ki\f{Xi, A2) given by 

K.ix.M ^ / I/ n(i + ^^)»-"-"'-"-""'^'-^". 

from which all the n-point correlation functions may be obtained by the usual deter- 
minant formulae of a matrix whose elements are the Kn{\i Aj). Defining the scaling 
variable y = N{Xi — A2), and using, in the large N limit, the expression of ( P38| ) for 
the product n^i(l + iy(^u~a ) )' ^^^r integration over t, one obtains 

-Im i — : — — ^ -e""^. (247) 



dXi IT I 27ri u + Gq{u) — Xi + ie 
Again one defines the pole u, and obtain 

dXi vr dXi 



^ ^ -Im(e-^"(^i-'^)) (248) 



ny dXi 
in which 

u(Xi - ie) = Ai - ReG'(Ai) - mp(Ai). (249) 
Therefore one ends up with 

i^^(Ai, A2) = -^e-^[^-^°^(^^)l sin[7ryp(Ai)]. (250) 
ny 
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Putting this expression of the kernel into the correlation functions, one sees that their 
only dependence in the eigenvalues of the external source is a scale factor through 
the density of state ; for intsnace the connected tow-point function is simply 



Pc(Ai, A2) = sin^[7rp(^l±^)y]. (251) 



We now consider the universality in the GOE and GSE ensembles in the presence 
of an external source. In the GOE case, we used a generalized HIZ integral which 
had two parts. One part was the t-integration, which gave simply 

* ^ ' - ^fa-J)(A.-A.) P^^' 
This leads now to the integral 

j{{[a,-it,)[a,-it,)j—^^e — — (253) 

In the large N limit, the saddle point equations for the tj{j = 1, 2) are 

t,- - iAi - ^7 y = (254) 



As in (|242|) , there is a saddle point ij, which behaves as ij ~ iXj in the large Xj 



domain. Then tj becomes 

ij = iXj - iReGiXj) + np{Xj) (255) 

Thus, this t-integration in the large N limit behaves as in the sourceless case, and 
d sin X 

yields — ( ) where x is iV7rp[(Ai + A2)/2](Ai — A2). More genrally, this part may 

be written as 

I=^K{x,y) (256) 

where K{x, y) is a kernel. This is due to the HIZ integral for (3 = 4. In our case, this 
kernel is a sine-kernel. 

For the u-integration, one has a product 

N 4 1 

n n -— - = exp[- E E - a,)] (257) 
7 j=i "i "7 7 j 
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There is a saddle point Wj, which behaves as ui = U2 — Xi and -us = ^4 ~ A2. Using 
these saddle points, the factor coming from this u-integral is the same as in (|179|) . If 



we denote this u-integral by I, it is easy to find that the derivative of I with respect 
to X is the sine kernel. More generally, we have 

-^I = K{x,y) (258) 

This is due to HIZ formula for f3 = 1 which gives 

E ^' 1 (259) 



In the sine kernel case, we have 

where the successive cancellation occurs in the higher order. 

Thus we have a universal two point correlation function in the presense of the 
external source, 

2,,-|- ,sinx,2 (i sinx, r°° smz 



p(Ai, A2) = p\m - C—^y - ;^(^) / ^dz] (261) 

where x = 7rA^p(A)(Ai — A2) and the external source appears only in the density 
of state p. The GSE case may be analyzed as the GOE and yields also universal 
correlation function with respect to the external source eigenvalues, as in ( p28|) . 



12 Universalities at the edges in GOE 

Near edges of the support of the density of states, it is well-known that a new scaling 



behavior takes place ( [pTj |29| . In the simplest case of the Wigner semi-circle the 
behavior is governed by an Airy kernel. Let us consider now the edge behavior of the 
two point correlation of characteristic polynomials, 

Fjv(Ai, A2) =< det(Ai - X)det(A2 - X) > (262) 
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when X is a real symmetric random matrix. If Ai and A2 are within the bulk of the 
support of the asymptotic density of states, we have found in the Dyson scaling limit 
that 

„ , , 1 d sinx , , 

i^^Ai,A2 =-1 263 

X ax X 

As we have discussed earlier , by using Grassmann variables, one finds for arbitrary 
Ai, A2 and finite N, 

FNiXi, A2) = / (det5)^e-™'+^^*'-^^ 

(tit2)^e-^*?-^*i-^^^*^^-2^^*^^n7^^^ - '!"h dtidh (264) 

(Ai - A2)"' (Ai-A2)'^ 

In the large N limit, the saddle points for the t^s are 



-iXa ± a/2 - A2 
ta = (265) 

The critical point corresponds to a degenerate quadratic form of fluctuations near 
the saddle point ; expanding then to next order (since there is a flat direction), one 
finds a new scaling limit when the parameters a and /3 are a = | and /3 = |, when 
the A's are at distance A^~" of the end point v^, and the taS at distance A^~^ of 
—i/^/2. Performing the scaling change of variables 

A, = V2-N-''xa 

ta = -^ + ^"^^- (266) 
the integral becomes, in this regime of large N, 

F«(A„A2) = /e-^-(-^^+-^^)+^-(-^--+----)[ ,^[^ ~ + ~ liM ridn 
J [Xi-M)^ (Ai-A2)^ 

= TT —-^[Ai"{xi)Ai{x2) - 2Ai'{xi)Ai'{x2) + Ai{xi)Ai"{x2)] 

[Ai — A2) 

(267) 

in which use has been made of the differential equation satisfied by the Airy function, 
Ai"(x) = xAi{x). Noting that 

(7^ T—)[Ai'{xi)Ai{x2) - Ai'{x2)Ai{xi)] 

0x1 0x2 

= -[Ai"{xi)Ai{x2) - 2Ai!{xi)Ai!{x2) + Ai{xi)Ai!\x2)] (268) 
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the second term of ( p64| ) is then simply the Airy kernel divided by (Ai — A2). One 
ends up with 



Fj^{Xi, A2) 




{xi)Ai{x2) — Ai'{x2)Ai{xi) 



Xi — X2 



(269) 



This confirms the result (|256|) , which was stated for more general kernels. 

Extending this analysis to ratios of characteristic polynomials, one obtains simi- 
larly an other term, expressible as an integral of the kernel. Thus we have obtained 
the two point correlation function near the edge for the GOE ensemble. The same 
argument can easily be transposed to the GSE. 

Let us now consider other edge problems that one meets when an external source 
matrix is added to the probability distribution. For instance one can tune the external 
source to create a gap in the spectral density of states p(A). At the critical point at 
which this gap closes, a new universality class appears and we have studied earlier 
the new scaling behavior at the origin p2|, 

It was found that the kernel Ki\[{x,y) for the CUE ensemble in the appropriate 
scaling limit, was 



For the GOE case, in the critical domain of the gap closing point, one finds a two 
point correlation function. 



13 Level spacing distribution in GOE and GSE 

The level spacing probability function E{s), the probability that there is no eigenvalue 
inside the interval [— s/2, s/2], is given by the Fredholm determinant 



K{x,y) 



^'{x)ij'{y) - 4>"{x)ij{y) - 4>{x)r{y) 



(270) 



x-y 




(271) 



det[l-K] 




(272) 
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Let us briefly review the derivation of tliis formula for E{s) This level spacing 

function is expressed by the Hamiltonian formalism. The derivation is due to Tracy 
and Widom p8| . 



In the GUE case, we have 

E{a,b) = det[l - K] (273) 

where we choose the interval (a, b) in general, and later we put (a, b) = (— |, |). The 
kernel K{x, y) acts on this interval, and we define 

K{x, y) = K{x, y)e{y - a)e{b - y) (274) 

in which 9{x) is the Heaviside function. By definition, the derivative of E with respect 
to the end points is 

^^^^ = -/?(M) (275) 

where 

k 

K = ^ (276) 

l-K 

,K is the Fredholm resolvant. One has a similar derivative with respect to a. Then, 
one obtains 

dEis) l.d 5 , , , ~,s s. / _x 

This leads to 

E{s)=eM- f ^(|'|)^^]- (278) 
where we write H{x,x) = K{x,x). For GUE, the kernel is the sine kernel, 

^^^^y^^mm^imvi. (279) 

x-y 

Operating with K on 0, one obtains 

q{x) =< x\-^\(P >,p{x) =< (P'\-^\x > . (280) 
L — K L — K 



This leads to 



K^^^y^ = ^(^)p(y)-^(y)p(^\ (281) 

x-y 
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Taking the derivatives of q{x) and p{x) with respect to b , one obtains a set of 
equations, which are a Hamiltonian system, 

Q = P(l-^),P = g(-^-l) (282) 

where Q{b) = q{h, —b : b) and P{b) = p{b, —b:b), q{b, —b, : b) , obtained by setting 
a = —b,x = b in q{b,a : x), (We have written exphcitly the interval dependence of 
q{x)). The Hamiltonian H which governs this dynamical system is simply 

H{b,b)=P' + Q'-^^ (283) 

b 

(P and Q have Hamiltonian form, Q = f^? -P = ~|^) with H = K{b, b)). For small 
s, we obtain from these two equations, 

7 o 65 o , . 

P = l + s + -s' + -s' + --- (284) 

3 

Q = -- — + ... (285) 
^ 2 48 ^ ' 

This leads to 

E{s) = 1 - s + 0{s^) (286) 

which agrees with the well known result 0. 

For the GOE ensemble, the kernel is a quaternion matrix, and the correlation 
functions are expressed by a quaternion determinant. 

The matrix kernel of GOE, a{x,y), is 

where we denote the sine kernel as s{x, y) = sm{x — y)/{x — y), and 

Js{x) = r s{y)dy - e{x) = £s{x) - e (288) 



where e{x) = ^sgnx. 

Tracy and Widom have found the Fredholm determinant for this matrix kernel 



28|. Their result for the interval (a, &) reads 

E{s) = exp[-l r H{z, z)dz T H{z, -z)dz] (289) 
2 Jo 2 Jo 
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for GOE. The Hamiltonian H is same as the Hamihonian for the GUE given in ( p83| ). 
We have 

= (290) 



For small s, one obtains from ( p89| ) 



E(s) = l-s + 4^' + 0(s^) (291) 

which agrees with the known result (We have dropped a factor vr in the sine kernel. 
The correct coefficient is 
For GSE, we have 

E(|) = 1 [exp[-i I\h{z, z) + H{z, -z))dz] + exp[-i j\H{z, z) - H{z, -z))dz]] 

(292) 

For small s, it gives 

E{s) = 1- s + 0{s^) (293) 

When the external source matrix A has only two dsitinct eigenvalues -a and +a, 
each of them y times degenerate, a gap in the spectrum around the origin may be 
created by tuning appropriately the parameter a, . At some critical value of a the gap 
closes and its vicinity leads to a new interesting universality class. For the GUE case, 
we have given in an earlier work |2^ the equation satisfied by i?(s)in the oppropriate 



scaling vicinity of the gap closing point. In this gap closing case, the problem is 
again governed by a Hamiltonian system with now three different Qi,Pi{i = 0, 1, 2), 
as shown in ||23|. They satisfy 



0. = H (294) 



The Hamiltonian H = K{b, b) reads 



H = bP2Qo + Q2Pi + QiPo-uPiQo-vP2Qi 

+ -^[Pi{b)Q,ia) - Q2{b)P2{a) - Qo{b)P^{a)] 
b — a 

X [Pi{a)Q,{b) - Q2{a)P2{b) - Qo{a)Po{b)] (296) 
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for the interval (a, b). 
The kernel K{b, —h) is 

^(^^ ^ giWPi(-6) - Q2{h)P,{-h) - Qo{b)Po{-b) ^297) 

Using the expressions of ( p^89| ) and ( [^92| ), we obtain the level spacing function both 

for GOE and GSE near the gap closing point. The second derivative of E{s) is the 

level spacing probability p{s), which is the probability density that two successive 

eigenvalues lie at distance s. For s small, it is easy to verify that p{s) is linear in 

s, a characteristic behavior of GOE. The small s expansion may be computed from 

1 

the expansions of Qn,Pn, which are Qo{b) = — r(-) + , & + 0{b'^),Qi{b) = 

An 4 zvT'^/^ 

-^nl)b + 0{b'),Q,{b) = -^ril) + 0{b),Po{b) = -^b' + 0{b'),Pi{b) = 
^ +0ib),P2ib) = ^b + Oib^). 



Thus we have obtained new results for the level spacing probablity of GOE. This 
may be important, for instance, in the discussions of universality for the the energy 



spectrum of quantum dots with interactions |2T|. Our study of the external source 
problem may be related to the questions of the distribution of cycles in the permu- 
tations with external source in the GOE case, and to the crystal growth in a 



random environment [19| or to the spin glass problem |20 



14 General p 

We have discussed hereabove the HIZ integration for arbitrary (3. The two point 
correlations of the characteristic polynomials are given for the general (3 case by 



< det(Ai-X)det(A2-X) >= C J f[{Xi-x,) l[{X2-x j) Y[\xi-x jlf^e-^^""' f[dx, 

i=l j=l i<j i=l 

(298) 



57 



in which the eigenvalues of X are Xi. In the (3 = 1 GOE case, this correlation function 
is simply 

^(Ai,A,) = = (299) 

V Zx'^ 2 X ax X 

In the GUE case, this two point correlation function becomes the kernel itself 



K]\f(Xi, A2), and it is given by W — Ji(x) 

V 2x 2 X 

It is easy to derive the general expressions for arbitrary (3 from the relation ( p.63| ). 
If we use the dual /?, which is equal to ^ and k = 1 in ( |163|) ???. 

The HIZ integral has a finite number terms when | — | is a half integer ; for 
instance for (3 = 1, the GOE ensemble is dual to the GSE, for which we know 
that the expansion is finite . From the saddle-point large N analysis, this two point 
correlation function of characteristic polynomials in the Dyson limit is 

""^Jl i(x) (300) 

The Bessel functions of half-integer order is expressed by 

^ ^^'^H-^n^y (301) 

For instance, when (3 = 4, k = 1, (GSE), we have shown that indeed it is a Bessel 
function Jo(x) in ( |142| ). 

The expression in terms of Bessel functions has been derived for the generalized 
Selberg integrals by Aomoto [0, and further studied by Forrester ^T[. It may be 



interesting to note the following fact. 

15 Summary 

In this paper, we have studied the correlations of the characteristic polynomials of 
random matrices, which are either real symmetrix (GOE) or quaternionic self-dual 
(GSE). It was shown that they are universal in the Dyson's limit with respect to 
an external matrix source linearly coupled to the random matrix. As usual the 
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correlation functions are only sensitive to the external source through a scale set by 
the mean spacing. 

For the ratio of characteristic polynomials, we have applied supersymmetric tech- 
niques, and obtained as a by-product the resolvent correlation functions of the GOE 
and GSE ensembles. The method required using group integrals which are no longer 
of the semi-classical type studied by Harish-Chandra, Itzykson and Zuber. However 
we have found that, in the limit required by Dyson scaling, these integrals can be per- 
formed both for the GOE and the GSE. We have also used an alternative method, in 
which one computes ratios of characteristic polynomials, and we have obtained again 
the correlation functions by a replica method, in the zero-replica limit. 

The level spacing probability has been studied for the GOE and the GSE ensemble, 
and we have given an explicit representation for E(s) in the closing gap case, obtained 
by tuning of the external source. 
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